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II.  SPATIAL  FEATURES  OF  ELECTRON  AND  PROTON  FLUXES  IN  THE 

NORTHERN  AURORAL  ZONE 


II.  A.  SPATIAL  FEATURES 

Electron  (0.19  -  3.2  MeV)  and  proton  (0.64  -  35  MeV)  data 
in  the  magnetic  media  provided  by  the  Japanese  collaborators 
appear  in  both  four  and  sixteen  seconds  intervals.  Using 
software,  all  the  data  -  counting  rates,  latitude,  longitude, 
altitude,  magnetic  field,  McLwain's  parameter  L,  local  time, 
magnetic  time,  and  the  angle  X  between  the  telescope  axis  and 
the  local  magnetic  field  -  were  averaged  to  find  the  per  second 
values  of  these  quantities.  Later,  average  of  these  quantities 
were  further  averaged  over  1-deg  latitude  bin.  In  the  next  step, 
satellite  revolution  numbers,  were  divided,  arbitrarily,  into 
satellite  pass  numbers  over  the  different  global  zones.  Further 
analysis  continued  from  this  point,  where  in  the  refined  data, 
the  counting  rates  are  in  per  second  per  degree  latitude. 


II.  A.  1.  SURVEY  PLOTS 

The  next  step  in  the  data  analysis  was  to  generate  plotts 
of  the  satellite  pisses.  In  each  plot,  the  average  counting 
rates  were  plotted  against  latitude  bins,  longitude,  altitude, 
magnetic  field,  L,  local  time,  magnetic  time,  and  the  angle  x> 
Fig.  1  is  such  a  plot.  These  plots  are  useful  to  locate  the 
position  of  the  peak  counting  rates,  and  the  peaks  in  different 
global  zones.  The  dotted  lines  indicate  protons  and  the  solid 
lines  electrons.  Zero  values  along  the  horizontal  axis  are  to  be 
ignored.  They  merely  denote  absence  of  data  following  the  end  of 
the  satellite  passes.  Relevant  literature  reports  of  the  peak 
positions  of  the  proton  aurora  at  L-4,  and  electron  aurora  at 
L<-5.  The  plots  seemed  to  show  peak  values  occuring  in  a  number 
of  L  values.  To  separate  the  auroral  zone  data  for  further 
analysis,  a  threshold  L  value  was  used.  From  a  careful  overview 
of  the  entire  data  set,  any  peaks  appearing  beyond  L-S  was  taken 
to  belong  to  the  auroral  zone.  Fig.  2  shows  the  B-L  vales  of  the 
peak  rates  falling  to  different  global  zones  in  the  northern 
hemisphere.  From  left  to  right,  the  first  cluster  of  points  are 
correspondii^  to  locations  of  the  peak  flux  in  the  equatorial 
zone,  whTch  is  then  followed  by  the  low-latitude  zone,  the 
mid-latitude  zone,  and  the  auroral  zone. 

II.  A.  2.  GLOBAL  PROFILE 

The  latitude  and  longitude  of  the  location  of  the  peak 
counting  rates  were  plotted  in  geomagnetic  coordinates  to 
generate  Fig.  3.  It  shows  the  global  profile  of  the  northern 
auroral  zone.  Every  square  represents  the  location  of  a  peak 
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Fig.  la  and  lb:  Identified  from  the  left  are  the  low-latitude,  the 
mid-latitude,  and  the  auroral  zones  of  particle  precipitation  in  the 
northern  hemisphere.  Dotted  curves  are  for  protons,  and  solid  ones  f 
electrons . 
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Fig.  2:  B— L  coordinates  of  the  peak  flux  profile  of  the  equatorial  zone, 
the  low-latitude  zone,  the  raid-latitude  zone,  and  the  auroral  zone  (from 
left  to  right) 


o  CP  o  cn  a  a 
Q  a  QD  OBD  OG  n  Q 
□  a  anna  oh  on  piam 
o  «Da  no  Hi  OBD  a  a  a  \ 
GO  an  am  0  0  oiisa  aoDoa 

Q  CO  CBCD  SI  ■■  HOOD 

■  cm  CX3  oconso 

a  a  □■■■■ 

n  a  a  □  a  ooms 

acs  a  oaao 

n  a  a  I 


□  a  OQ  Qo 

0  a  mo  o  D 

a  ODO  mai 

a  DO  mm  o  on 

na  o  ooo  o  ooioksd  mdi 

■m  m  am  amaD  oasaBa 


0.0  36.0  72.0  108.0  144.0  180.0216.0252.0288.0324.0360.0 

GEOMAGNETIC  LONGITUDE  (DEG) 

Fig.  3.  Global  peak  flux  profile  of  the  northern  auroral  zone 
in  geomagnetic  coordinates. 
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counting  rate.  The  global  profile  runs  parallel  to  the  minimum 
magnetic  field  profile. 


II.  A.  3.  SCATTER  PLOTS 

Scatter  plots  are  very  important  to  study  the  dependences 
of  the  peak  counting  rates  on  one  or  more  of  these  quantities  - 
longitude,  altitude,  L,  local  time,  x  angle,  etc.  Software  was 
developed  to  make  scatter  plots  of  the  peak  rates  vs  the  above 
quantities,  and  also  one  or  more  of  the  above  quantities  vs  the 
rest  of  the  quantities.  These  have  been  used,  particularly,  to 
investigate  the  dependences  of  the  peak  rates  upon  more  than  one 
variables.  Fig.  4  shows  a  scatter  plot  of  the  peak  rates  of 
protons  vs  the  x  angles. 


II.  A.  4.  INSTRUMENT  RESPONSE  FUNCTION 


The  S-1  telescope  response  function  varies  with  the  x 
angle.  Fig.  5a  through  5f  shows  the  response  function  for 
different  orientation  of  the  telescope  axis  with  respect  to  the 
local  magnetic  field  direction.  The  pitch  angles  (angle  between 
the  velocity  vector'  of  the  particles  and  the  magnetic  field 
direction)  of  particles  are  given  along  the  abscissa  of  the 
plots.  The  response  function  is  maximum  for  particles  of  pitch 
angle  a-90"  when  x*90^.  Gradually  decreasing  X  has  gradually 
decreasing  response  function.  What  is  worthnoting  is  the 
correlation  between  x  and  response  function  f. 


II.  A.  5.  SUPERPOSED  PLOTS 

Software  was  developed  to  superpose  satellite  passes  whose 
peaks  fall  in  a  certain  longitude,  altitude,  L,  local  time,  and 
beta  angle  bins,  for  in-depth  study  of  the  "weather  and  climate" 
of  the  auroral  zone.  The  software  developed  makes  peak-to-peak 
superposition  of  the  satellite  passes.  The  off-peak  rates 
falling  in  a  particular  latitude  bin,  are  averaged  over  the 
number  of  passes  contributing  non-zero  rates  in  that  bin.  The 
plots  of  superposed  passes  show  a  larger  number  of  latitude  bins 
than  the  actual  range  of  latitude  for  the  location  of  the  zone. 
This  is  due  to  the  fact  that  different  passes  have  peaks  at 
different  latitudes  in  the  auroral  zone.  Latitude  bin  marked  by 
50  is  the  position  of  the  peak  rates.  Each  latitude  bin  is  1" 
wide . 


RESPONSE  EUNCTION  RESPONSE  FUNCTION 
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II.  B.  1.  LONGITUDE  DEPENDENCE 


Longitude  range  was  binned  in  36°  intervals.  Satellite 
passes  having  peak  positions  in  a  longitude  bin  were  superposed 
peak-to-peak .  Fig.  6a  is  such  a  plot  in  the  longitude  range  0°  - 
36°  for  protons.  Fig.  6b  shows  the  longitude  variation  of 
protons  (dotted  line)  and  electron  (solid  line)  fluxes.  Minimum 
flux  occurs  in  the  longitude  range  36°  -  72°.  The  flux  then 
increases  and  reaches  maximum  between  216°  -  252°.  From  the 
scatter  plot  of  x  vs  geomagnetic  longitude  in  Fig.  6c,  it  is 
seen  that  the  maximum  response  function  at  the  longitude  range 
of  the  maximum  counting  rate  is  not  the  factor  to  enhance  the 
counting  rate.  In  the  said  longitude  range,  the  instrument 
response  function  uniformly  covered  the  response  function 
corresponding  to  X“  90+35°,  instead  of  X"  90°  alone.  Thus 
auroral  proton  flux  does  show  longitude  dependence. 


II :  B.  2.  ALTITUDE  DEPENDENCE 


Altitude  range  (350  -  850  km)  was  binned  in  intervals  of 
100  km.  Peak  counting  rates  falling  in  a  specific  bin  were 
superposed  peak-to-peak.  Fig.  7a  is  such  a  plot.  Later  the  peak 
fluxes  have  been  used  to  study  the  altitude  dependences.  Fig.  7b 
shows  the  altitude  dependences  of  both  electrons  (solid  lines) 
and  protons  (dotted  lines).  Electron  flux  show  a  strong 
exponential  dependences.  Proton  flux  also  show  some  altitude 
dependences  up  to  750  km.  The  exponential  increase  of  electron 
flux  cannot  be  due  to  longitude  dependence  or  the  maximum 
response  function  of  the  instrument. 


II.  B.  3.  DEPENDENCES  ON  X 


Because  of  the  correlation  between  x  and  the  response 
function,  it  is  expected  that  a  plot  of  x  and  the  counting  rates 
should  show  the  maximum  counting  rate  at  X"90°  when  response 
function  is  the  maximum,  and  gradually  diminishing  counting 
rates  as  X  decreases  that  is  as  the  response  function  decreases. 
To  investigate  the  expected  variation  of  particle  flux,  the 
range  of  X  was  binned  in  10°  intervals  and  the  superposition  of 
passes  was  made  in  each  of  the  bins.  Fig.  8a  shows  one  of  the 
plots  of  superposed  passes.  A  stray  counting  rate  peaks  at 
latitude  bin  70.  The  peak  counting  rates  from  such  plots  were 
used  to  plot  Fig.  8b  for  protons,  and  Fig.  8c  for  electrons.  It 
is  seen  that  the  flux  variation  follow  sine  function  of  X*  This 
is  probably  due  to  the  fact  that  the  pitch  angle  distributions 
of  auroral  zone  protons  and  electrons  can  be  described  by  some 
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Fig.  6b.  Illustrates  the  longitude  dependence  of  the  peak  proton 
(dotted)  and  the  peak  electron  (solid)  fluxes. 
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Fig.  7a.  Plot  of  superposed  passes  having  peaks  in  the  altitude 
range  750  -  850  km. 


Fig.  7b.  Illustrates  altitude  dependence  of  both  electron  (solid) 
and  proton  fluxes  (dotted). 
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Fig.  8a.  A  plot  of  superposed  passes  having  peak  positions  within  Xi= 
70  -  80  degrees. 
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Figs.  8b  and  8c.  Illustrates  the  variations  of  the  peak  flux  as 
a  function  of  Xi.  Fig.  8b  (top)  for  protons,  and  Fig.  8c  (bottom) 
for  electrons. 
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powers  of  sine  function. 


II.  B.  4.  DEPENDENCES  ON  L 


Contrary  to  the  reported  information  available  in  the 
literature,  the  variations  of  the  peak  fluxes  of  electrons  and 
protons  in  L  space  depend  on  the  pitch  angle  of  the  particles, 
as  are  evident  from  the  plots  of  the  peak  flux  as  a  function  of 
L  for  different  X  bins.  These  are  illustrated  by  Figs.  9a 
through  9j.  Both  electron  and  proton  (L  <  6)  fluxes  show 
exponential  dependences  on  L.  Proton  flux  shows  a  second  peak  at 
L-6.  At  lower  X  (<  60*)  values  which  is  true  for  lower  pitch 
angles  of  the  particles,  this  behavior  changes.  Exponential 
dependence  on  L,  possibly,  contributes  to  the  exponential 
dependence  on  altitude  h. 


II.  B.  5.  WIDTH  OF  THE  AURORAL  ZONE 


All  the  satellite  passes  for  proton  and  electron  counting 
rates  were  superposed.  Fig.  10a  shows  the  plot  of  the  superposed 
passes  for  electrons  and  Fig.  10b  that  for  protons.  It  is  seen 
that  the  proton  flux  zone  is  wider  than  the  electron  flux  zone. 
The  Full-Width-at-Half-Maximum  (FWHM)  for  the  electron  zone  is 
~  5®  and  that  of  the  proton  zone  is  ~  7". 


III.  TEMPORAL  FEATURES  OF  ELECTRON  AND  PROTON  FLUXES  IN 
THE  NORTHERN  AURORAL  ZONE 

Both  the  electron  and  proton  fluxes  show  short  time  and 
long  time  variations. 


III.  A.  DEPENDENCES  ON  LOCAL  TIME 


The  local  time  was  binned  into  eight  intervals  of  three 
hours  eaph-  As  before,  satellite  passes  having  peak  rates  in  a 
particular  bin,  was  superposed  peak-to-peak .  Fig.  11a  is  one  of 
the  plots  of  the  superposed  passes.  The  peak  fluxes  in  the  time 
bins  were  used  to  investigate  the  local  time  variation 
illustrated  in  Fig.  lib.  The  flux  is  seen  to  peak  to  different 
values  at  two  different  local  times.  A  scatter  plot  of  local 
time  vs  X  shows  that  these  peaks  are  due,  partly,  to  the  maximum 
instrumental  response  functions  corresponding  to  X  ~  90*. 
However,  since  the  early  morning  hour  peak  is  higher  than  the 
afternoon/evening  peak,  we  conclude  that  the  flux  becomes 
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Fig.  9b 
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maximum  at  the  early  morning  hour.  A  similar  variation  is 
observed  for  magnetic  local  time. 


III.  B.  LONG-TERM  VARIATION 


Monthly  average  of  the  electron  and  proton  fluxes  were 
made  by  the  peak-to-peak  superpositon  method.  Fig.  12a  shows  one 
of  the  monthly  averages.  Plot  of  these  monthly  averages  shows 
periodic  variation  of  electron  and  proton  fluxes  -  one  of  them 
become  relatively  much  higher  than  the  other  one,  as  illustrated 
in  Fig.  12b.  At  times  (viz.,  October,  1984),  electron  fluxe  can 
be  as  much  as  9  times  that  of  proton  flux.  It  is,  also,  seen 
that  these  fluxes  show  a  kind  of  periodicity  over  a  number  of 
months.  The  first  maximum  occurs  in  May,  1984,  then  in  October, 
followed  by  February  and  August  of  1985,  and  then  in  January 
1986.  Monthly  mean  of  geomagnetic  storm-time  index  Dst  (Fig. 
12c),  does  not  show  any  correlation  with  the  monthly  flux 
variation.  It  has  been  found,  at  the  preliminary  investigation, 
that  the  flux  enhancement  is  not  entirely,  due  to,  response 
function,  L,  altitude,  longitude,  or  local  time  variations.  This 
and,  at  least,  the  relative  variation  may  be  due  to 
magnetospheric  conditions.  The  investigation  is  in  progress. 


IV.  PUBLICATIONS/PRESBNTATIONS  HADE 


One  article  has  been  published  in  the  Journal  of 
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1295-1306,  1993.  The  title  is  "Solar-induced  variation  of  proton 
precipitation  near  the  equator".  This  article  will  also  appear 
in  abbridged  format  in  the  proceedings  of  the  Solar-Terrestrial 
Physics  Workshop  held  in  Canada  in  1992.  Another  article  is  in 
press  for  publication  in  Advances  in  Space  Research  of  COSPAR 
journal.  The  title  is  "Significant  variation  of  proton 
precipitation  near  the  equator".  This  article  was  presented  in 
the  World  Space  Congress  Meeting  held  in  Washington,  DC,  in 
August,  1992. 


V.  PROJECTED  PUBLICATIONS/PRESENTATIONS  FOR  1993-94 


Two  articles  are  almost  ready  for  submiting  to  journals 
for  publications.  One  of  them  is  on  the  midlatitude  zone,  and 
the  other  on  the  auroral  zone.  A  paper  on  the  particle 
precipitation  in  off-equatorial  global  zones  has  been  submitted 
for  presentation  in  the  7th  Scientific  Assembly  of  the 
International  Association  of  Geomagnetism  and  Aeronomy  to  be 
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Fig.  12a.  Plot  of  superposed  passes  having  peaks  in  the  month  range  numbered 
8  and  9  (i.  e.,  for  August)  of  1986. 
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Fig.  12b.  MontRly  variations  of  peak  electron  (solid)  and  proton 
(dotted)  fluxes. 


MONTHS 

Fig.  12c.  Monthly  mean  values  of  Dst  index  for  the  years  1984-86. 
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held  in  Buenos  Aires,  Argentina,  in  August  8-20,  1993.  Another 
article  will  be  presented  in  the  Western  Pacific  Geophysics 
meeting  to  be  held  in  Hong  Kong  in  July,  1994. 


V.  PROJECTED  TASKS  FOR  1993-94/REQUEST  FOR  RENEWAL 


Investigation  of  the  global  zones  lying  to  the  south  of 
the  geomagnetic  equator  (lower  part  of  Fig.  13)  has  not  been 
done  yet.  Studying  the  spatial  and  temporal  features  of  the 
particles  precipitating  in  these  zones  will  be  exciting. 
Interesting  features  lie  in  the  comparative  study  of  the 
southern  and  northern  zones.  Software  prepared  to  study  the 
northern  global  zones  will  be  useful  with  slight  modification  to 
study  the  southern  global  zones.  However,  these  regions, 
specially  the  one  falling  in  the  low-latitude  zone,  is  expected 
to  be  complicated  by  the  South  Atlantic  Anomaly.  About  five 
man-months  are  expected  to  complete  the  work. 

I  am  making  the  second  and  the  last  renewal  request  in  the 
interest  of  uncovering  facts  from  southern  global  zones. 


Fig.  13.  Illustrates  the  global  peak  flux  profiles  of  the 
equatorial  zone,  the  northern  (upper)  low-,  mid-,  and  auroral 
zones,  and  the  southern  (lower)  low-,  mid-,  and  auroral  zones 


